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SUMMARY 


The investigation presented herein of the sound field about a point 
source over a plane boimdary in the presence of a vertical temperature 
gradient represents one phase of a general program of research in atmos- 
pheric acoustics. 

In this report geometrical ray acoustics is employed to derive the 
ray paths and the intensity distribution about a source located in a 
viniformly stratified mediimi. A comparison of the relative effect of the 
stratification on the intensity from a directive source and from a spher- 
ical source is made. 

Also presented is a theoretical analysis of the sound field in the 
"shadow zone" (diffraction region) formed over an absorbing plane boundary 
in a temperatiire -stratified atmosphere. The analysis holds for both two 
and three dimensions. The boiondary condition at the plane is given by a 
normal acoustic impedance independent of the angle of incidence. As in 
the corresponding problem of underwater sound where the boundary is a 
pressure -release surface, it is found that the major portion of the soimd 
intensity in the shadow region decays exponentially with distance at a 
rate proportional to the one -third power of frequency and to the two- 
thirds power of temperature gradient. The effect of ground impedance 
enters mainly throiagh its resistive component. The rate of sound decay 
for a pressure-release boundary (zero impedance) is found to be 2.5 times 
that for a reflecting boundary (infinite impedance). 

Finally the results of measiirements of the sound distribution in a 
two-dimensional laboratory chamber in which a large temperature gradient 
was created are reported. Measurements were made at various frequencies 
for both reflecting and absorbing boundaries, and the results aire found 
to be in good agreement with theory. 
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INTRODUCTION 


A study of the propagation of soimd in the atmosphere must take 
account of the fact that the medium is^ in general, neither homogeneous 
nor quiescent. The principal atmospheric factors affecting soiind propa- 
gation seem to he humidity, turbulence, and wind and temperature gra- 
dients (ref. 1). In addition to these, the effect of the gro\md must be' 
considered (ref. 2 ). With the exception of humidity, the effect of which 
is generally small, these factors are clearly interrelated, the strength 
of one being somewhat dependent on the presence of the others. 

The presence of tiirbulence causes attenuation of sound as a result 
of a scattering of the sound energy, the attenuation increasing with fre- 
quency and wind velocity. Turbulence has generally been considered the 
factor of major importance among those mentioned (ref. 5), although no 
systematic study of it seems to have been made. The existence of wind 
and temperature gradients is evidenced in the continuous refraction of 
the acoustic rays or energy paths. This continuous refraction in itself, 
has little influence on the sound intensity; however, the effect of refrac 
tion in attenuating soimd may be greatly enhanced, often exceeding that of 
t\arbulence, owing to a peculiar influence of the groimd boundary resulting 
in "shadow formation" as illustrated in figiare 1 . This phenomenon, which 
has been rehlized and studied in some detail in the field of underwater 
sound, seems to have been overlooked in the atmosphere. However, unpub- 
lished field measurements obtained at the I'to.ssachusetts Institute of 
Technology have clearly shown the importance of shadow formation, and 
the present detailed study of the sound field in the temperature— created 
shadow zone was si;iggested by these observations. 

This investigation was conducted at M.I.T. under the sponsorship 
and with the financial assistance of the National Advisory Committee for 
Aeronautics . 


SYMBOLS 


b frequency parameter, PCq^'ToO) 

c sound velocity, function of z 

Cq sound velocity at ground 


H "height" of atmosphere 

Hn^^^(x) Hankel function of second kind of order n 
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h 

source height 

I 

intensity 

Jn(x) 

Bessel fiinction of order n 

k 

propagation constant 

L 

total power output 


- k2 

r 

radial coordinate, distance from source 


distance from source at height h to shadow boimdary at 
height z 

To 

absolute temperature at ground 


travel time from soiorce at height h to shadow boundary at 
height z 


V = Uoe"^” 


Z 

acoustic impedance of ground 

Zi/5(u) 

linear combination of one-third order Bessel functions 

z 

vertical coordinate, receiver height 

2o 

first-order zero of Q, Q(2 q) = 0 

a 

so;md velocity gradient 

P 

temperatxnre gradient 
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^ impedance ratio of gro\md, Z/pqCq 

T] admittance ratio, l/^ 

6 resistive component of impedance ratio 

p density of air 

Pq density of air at ground 

T defined "by k = ^ cos t 

0 velocity potential 


RAY ACOUSTICS 


For many problems involving the sound field about a source in a 
nonhomogeneous medium the high-frequency approximation afforded by geo- 
metrical (ray) acoustics is sufficient. Unless only qualitative results 
are desired, however, the method becomes impractical in the neighborhood 
of reflecting boimdaries, and it clearly fails inside a shadow zone. In 
the case of a sound field from a spherical soiirce, ray acoustics predicts 
very small effects from the temperature gradients which are normally 
enco\mtered in the atmosphere. Thus, for instance, the normal adiabatic 
lapse rate of 1° C per 100 meters causes a sound diminution in the hori- 
zontal direction of only 10“^r^ decibels (r in kilometers) in addition 
to the diminution indicated by the inverse -square law for short ranges. 
There is, however, one case amenable to ray acoustics where a ten 5 )erature 
gradient can have quite a significant influence on a soxmd field and that 
is when the source is directive. 


The path of a soixnd ray is determined directly by Fermat's principle 
which requires that the path be such that the time of traversal is a 
minimijm (or extremtim) . This condition is met by requiring that the inte- 


gral 



be stationary or, expressed symbolically, that 



0 , 


where ds is an element of the path, c is the velocity at that point, 
and 6 is the conventional symbol for the variation of a functional. It 
is assumed that c is a fimction of one coordinate only, namely, c(z), 

so that ds = ^1 + (r ' )^ dz, where r' = dr/dz. The Euler eqmtion corre- 
sponding to the problem 
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6 


Vl + (r')^ 

c(z) 


dz 


= 0 


is 


whence^ 



(r’)^ 

1 + (r')^ 


where A is a constant. The general equation for the ray path is then 


r 



dz 


c 


2 


( 1 ) 


Putting 


— = tan 0, where 
dr 


9 is the angle of inclination of the ray, gives 

A = shL 

COS 0 


which is Snell's law. Thus^ in the integrand in equation (l)^ A plays 
the role of a ray parameter. The integration of this equation yields 
circular paths for a constant velocity gradient and cycloids for a veloc- 
ity gradient proportional to the square root of one coordinate (constant 
temperature gradient). 


Distance to Shadow Zone 

Consider a source to he located a height h above ground (at z = O) 
and suppose the sound velocity to diminish with height. Those rays which 
just graze the earth, that is, which have zero slope at z = 0, then define 
a bounding surface (in three dimensions) separating the sound region from 
the shadow region. For these rays, clearly the ray parameter A 
equals c^, with c^ the velocity at the ground. The horizontal dis- 
tance between a source located at z = h and the start of the shadow zone 
is then 


h 


c(z) dz 

- c(z)2 


r(h) = 


0 


( 2 ) 
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If the receiver is located at a height z on the shadow boundary^ then 
the horizontal distance from source to receiver is 

Tq = r(h) + r(z) 

For a constant gradient of velocity, c(z) = Cq - az and 



For a constant temperature gradient 



(3) 



and 



The approximations hold for small values of the gradient. Equation (A) 
has been plotted in figure 2 for the case in which source and receiver 
are at the same height. 


The travel time tQ of the ray from source to shadow boundary is 
discussed subsequently. This time may be broken up into t^ = t(h) + t(z), 
where 



The value of r' has been substituted from eq\aation (2). For a con- 
stant temperature gradient, 



( 5 ) 


Intensity Distribution 

Once the paths of the sound rays have been determined, it is a simple 
matter, as shown in reference 5^ 'to calculate the intensity distribution 
by considering the geometrical divergence of these rays. Thus, consider 
a spherical source with a total power output L. The intensity on a 
small sphere of radius a centered on the source is then L/4«a^. Con- 
sider now two neighboring rays leaving the sovirce at angles of elevation 
0Q and 0Q + 600* Then, through the space bounded by the two surfaces 

fonned by rotating these rays about the z-axis, a certain amount of power 
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5L flows eqml to that which falls on the sphere zone cut out by these 
rays; that is. 


6L 


L 


2rta^cos 0 q60q 


L 

- cos 


0O&0O 


The intensity I at any point (r,0) is then the ratio of this power 6L 
to the area 2rtr6a thro\:igh which it flows at this distance, 6a being 
the perpendicular separation of the rays. Thus, 

I = ^ ^ (6) 

4jtr 6a 


The value of the interval 6a along the ray remains to be calculated. 
From figure 3 i't can be determined that 


6a = sin 06r 
= sin 06 / dr 


= sin 06 / cot 0 dz 


= -sin 


.Q n 


60 


p 

sin 0 


dz 


The q-uantities 50 and dz appearing in the integrand may he computed 
from Snell law, which reads 


cos 0 __ c(z) 

cos 0O Cq 

cos(0 + 50) __ c + 5c 

cos(0q + 60o) 

Neglecting 5c and expanding the cosine gives 

(tan 0)60 = (tan 0 o)50q (T) 

Differentiating Snell *s equation gives 


dz 


c d0 
c 'tan 0 


( 8 ) 
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dc 

where c' = — . Combining equations ( 7 ) Q-nd (8) gives finally 
dz 


6a = -sin 0 tan 0 q50q 


c d0 
c ^ sin^0 


(9) 


These results are now applied to the case of a constant temperature 
gradient. For simplicity the source and receiver are assumed to be at 

( ’ 1 02 

1 + -^1 y 

^0/ 


p 

/cos 0o\ 


Vcos 0 j 


and equation (8) becomes 

Tq sin (/) sin r ^ 


50 ^ _2 

P cos^0o ^ 


cof^0 d0 




cos 


500 L 


sin(0o - 0) - (0o - 0)sin 0 sin 0 q 


The value of this quantity at the so\irce height is obtained by putting 
0 = -0Q (note that the above integral may not be evaluated directly 
from 0Q to -0Q since it diverges); whence^ 


60 ^ sin 0 o 

^0o P cos^0Q 




( 0 = “ 00 ) 


( 10 ) 


The horizontal distance r from source to receiver (at the source height) 
may be expressed in terms of 0^ as 


r 


cot 0 dz = 




sin 200 ) 


( 11 ) 


Here, use has been made of the differentiated form of Snell's law, equa- 
tion (8). In terms of Iq, the intensity from the same source in a 

homogeneous mediiom ^ , the excess attenuation due to the inhomo- 

geneity is obtained from equations (6) to (lO) as 
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j _ 1 200 + sin 200 

Iq h tan 0Q^1 + 0otan 0o) 


( 12 ) 


where 0o is related to r in equation (ll) . For small values of 
3r/To eqiiations (ll) and (l2) combine to give 


lo 


1 


12 



2 


For a soiarce I' with a directivity pattern f 

(jS ^ 


(15) 


The relative effect of a temperature gradient on a spherical source and 
on a directive one is illustrated in figure k, where equations (l2) 
and (15) have been plotted as a function of the distance parameter pr/To 

from equation (ll) for a simple soiorce. The second curve is for a dipole 
source oriented at 45° to the horizontal, f“(0o) = - (it/4^ . It 

is seen that out to a distance of « 1.8 (about 5 kilometers for the 

io 

normal adiabatic lapse rate) the effect of the temperature gradient is to 
increase the sound received from the dipole source but to diminish that 
from the spherical one. Unless a temperature inversion prevails, a 
shadow zone may, of course, begin to form before this point. 


DIFFRACTION INTO SHADOW ZONE 


In order to derive an expression for the sound distribution within 
an acoustic shadow zone brought about by a temperature stratification of 
the medium above a plane boundary, several approaches are discussed and 
evaluated. In one approach the formal Green's function solution to the 
wave eqxiation in such a medium is derived as a series of normal modes , 

The fmdamental but involved method of contoiu integration presented by 
Pekeris is discussed briefly and compared with the Green's function 
method, the case of a medi\im with a constant (acoustic) velocity gradient 
being used as a basis for comparison. Finally, an approximate solution 
is presented for the case in which a constant gradient of temperature 
exists in the medium. The rapid convergence of this solution in the 
shadow zone is such that the behavior of the sound field within this 
region is adequately expressed by the first mode in the series. 
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Formal Solution 

Consider a spherical source of sound which is driven at a constant 
frequency f = 03/2 tc, located a distance h above a ground boundary at 
z = 0, and assumed acoustically described by a normal impedance Z, 
independent of the angle of incidence. Suppose the medium to be strati- 
fied so that the sound velocity is a function of the one coordinate z. 

The sound field may then be derived from a velocity potential 0e^^ 
satisfying the wave eqixation 


V20 + 0 = 0 

c(z) 


(14) 


outside the source. 

The use of this equation entails two approximations from the outset^ 
namely^ the usual one that small amplitudes are required and, also, that 
the fractional change in any property of the medium must be small in the 
distance of a wavelength (see appendix A). The latter condition is a 
restriction to high frequencies but is not stringent even for the very 
high temperature gradients which were created artificially in the model 
weather chamber to be described subsequently. For example, if 6X is 
written for the change in wavelength in the distance of a wavelength, 
then the requirement (5A/A) « 1 is equivalent to (a/f) « 1 where 

dc 

a = — is the velocity gradient. In practice a is of the order of 
dz 

0.01 second"^ for a stable (adiabatic) atmosphere while, in the chamber, 
values as high as 120 second"^ were produced. 

Equation (l4) can always be separated in cylindrical coordinates 
by putting 0(r,z) = R(r)F(z), which gives 



where is the separation constant. Now, if it is assumed that 

(l) the velocity c(z) is a continuous monotonic function of z such 
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that c(o) = Cq and c(h) = 0, that (2) at z = 0 the boundary condi- 
tion is — ^ and that (5) for large values of z the solution 
dP/dz izjop 

represents upgoing waves, a set of eigenfunctions Fjj^(z) and corre- 
sponding eigenvalues my he dfetermined based upon these assumptions. 

The solution to the problem may then be expressed in the following equa- 
tion as an infinite sum of these eigenfunctions with undetermined 
coefficients: 


0 

m 

Here, Ho^ 2) (kmr) is chosen as the solution to equation (15) in order 

to insure outgoing waves in conjunction with the time factor In 

order to evaluate the coefficients A^ equation (l4) must be rewritten, 
introducing a source term on the right corresponding to a point source 
at r = 0, z = h, as follows: 


v2jl + -a£_ 0,-2 s(r)s(z - h) (18) 

c(z)2 

Here, 5 stands for the "delta" function. Substituting equation (17) 
into equation (18) and making use of equation (16) lead to 



ii, (a) 

dr ° 


r — H„^ (l^r) 


= - I &(r)&(z - h) 


Now, multiply both sides of this equation by Fj^(z) and integrate over 

the surface of a small cylinder enclosing the sotirce. Because of the 
orthogonality of the fractions, it is possible to obtain the following 
equation: 


-knF^ih) 


^H 


Fm2(z) dz 


= 2rtAjnr ^ 
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As r— >0, 


-4rtFjjj(h) 



> 

dz 


The solution is, then. 


0e 


io)t 


-ijte 


icDt 


00 


m=l 


Ho 





z) 


(19) 


where 


<iF’m/dz 


— ^ at z = 0. 
ia>p 


In two dimensions 


(x,z) the solution follows the same pattern; thus. 


0 


-2rti 

m=l 


e-^^l''‘Fjh)Fjz) 

'-'0 


(20) 


with the same boundary conditions on F. 

* 

The functions depend on the velocity profile c(z) and the 

boundary conditions at the ground. Two cases are now considered, one 
in which the velocity c(z) diminishes proportionally with the height 
(constant velocity gradient) and the other in which it diminishes with 
the sqiiare root of the height. 


Medium With Constant Velocity Gradient 

Equation (l^) has been solved in reference 6 for the case of a 
source located underwater with c(z) = az. In accordance with the nota- 
tion of reference 6, reverse the sign of the z -coordinate and assume the 
source to be located at z = cT and the groimd at z = 7. The eigen- 
function solutions to equation (16) are then 


Fjz) = zl/2j.^(ikmz) 
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where 






1 + 2 i ^ 
a 


J'(xni) = 0 


( 21 ) 


PQ^O 


Here n . 
equation (l 9 )> then takes the form 


q = — - — , and Xm = ikjjjy. The general solution. 


0 = - 2 ijt 



(kmr 

) Jin( ) *^in( ) 

[1 + (n^/xm^) 

]jin^(Xm) + 1 

(Jin') 

l^(Xm) 


( 22 ) 


For a pressure release at the boundary, t] — ^ and Jin(^) “ 

which case equation (22) comes under Pekeris ' solution. 

Discussion of method of solution . - It is perhaps worthwhile at this 
point to give a brief simimary of the method used by Pekeris to arrive at 
his solution in order to compare it with the method of normal modes 
employed herein. Pekeris starts from the wave equation (without a sovurce 
term) and divides the medium into two regions, one above the soxurce and 
the other below. The function F is expressed as a linear combination 

of the fxjnctions z^Z^Ij^j^(kz) and z^Z^j^^(kz) in each region, I and 

K being the modified Bessel fxmctions. The foxu* axbitrary constants 
involved are then determined (l) by the boundary condition at the water 
surface, ( 2 ) by the requirement that the pressure be continuous across 
the interface of the two regions, ( 5 ) by the requirement that at great 
depths the solution represent a downgoing wave, and (4) by the source 
strength. In this way the solution is expressed as an integral over the 
separation constant k as 


pOQ 

0 = / Jo(kr)F(z)k dk (25) 

the path of integration being determined to satisfy the source condition, 

that is, such that - J' 0 dS = 4«q, where q is the source strength. 

An investigation of the behavior of the integrand in the complex k-plane 
reveals that it is free of branch points but that it possesses poles 
and that these are all situated in the fourth qxiadrant. The integral 
in equation ( 25 ) is then split into two integrals by writing 


and the path of integration of the 


Jo(^) = i 




+ H, 


( 2 ) 


(kr) 


new integrals is deformed to enclose the first and fourth quadrants, 
respectively. This procedure results in the solution being expressed 
as a sum of residues taken at the poles of the integrand in the fourth 
quadrant together with two integrals taken respectively over the posi- 
tive and negative imaginary axes. These two integrals can be shown to 
cancel each other, and the remaining sum of residues is equivalent to 
equation (22) with Jin(^m) =0 q = oo. 


This method of solution is clearly tedious since for each case it 
must be ascertained whether branch points exist, that the integrals taken 
along the infinite axes vanish, and that the remaining integrals over the 
imaginary axes cancel. From the theory of normal modes these conditions 
would be expected to hold automatically. That the theory does not pre- 
clude the possibility of the existence of brstnch points but only requires 
that they shall not contribute to the solution is illustrated by the case 

where c(z) = Cq^( 1 + az which is also discussed by Pekeris. The 

integrand corresponding to equation (23) now contains three branch 
points, and the solution includes an additional branch line integral. 
However, by a detailed examination of the phase changes of this branch 
line integrand as it is taken around the branch cut Pekeris is able to 
show that its contribution is zero. The use of Pekeris* method, then, 
would seem to be dictated only in cases where the normal -mode theory 
breaks down. This occurs whenever the function F or its derivative is 
discontinuous, for then the necessary orthogonality relation which nor- 
mally results from equation (16), namely, 

(Vn - Vn)^ " - 0 = / Vn (2>t) 

^ ' z=0 ^ ^0 


no longer holds. In practice the first derivative must, of course, 
always be continuous (continuity of velocity), but a discontinuity in 
the density will make F discontinuous (continuity of pressure). Again, 
if there is a discontinuity in the sound velocity c(z), the eigenfunc- 
tions Fjj^ do not form a complete set (ref. 7) although they are still 

orthogonal; however, equation (24) still holds since only the second 
derivative is affected. In all of these cases the solution must include 
a branch line integral. In a homogeneous (infinite) medium this will be 
the entire solution since then there are no modes (continuous spectrum). 

Asymptotic solution . - In an appendix to reference 6 Pekeris derives 
an asymptotic expression for the imaginary -order Bessel functions 
appearing in the solution of equation (22) herein which is valid for 
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high frequencies ^ ^ largej . Use is made of this expression to 

arrive at a general relation for the roots of equation (2l) valid for 
an arbitrary boimdary. 


In the limit of large values of n 
asympt ot i cally 


reference 6 shows that 



where 


0 


in It 



(25) 


T 


-1 


1 4 - 



Since the value of T is small (at the roots it is proportional to 

0 1 /"X. 

— Q presented subsequently and it varies as f ' for high fre- 
0) ° 

quencies), the inverse tangent may be expanded as follows: 


0 


in^i 


i + . 

5 


Retaining only the first term gives 

e « 1 inT^ = - i n|^ + 


or 

X = ik 7 = in + A^n^/^e^^/^ + . . . (26) 

mm' m 

where This expression holds only at the boundary % - y , 

since elsewhere T need not be small. 

Now, the Hankel fimction appearing in equation (22) is represented 
asymptotically (large value of r) by 


i 6 
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Since the roots from equation (26) have a positive imaginary part, this 
clearly leads to an exponential decay at large distances from the source. 

The rate of decay in decibels is ^8.68 ^ A^n^/^^y/7 decibels per meter 

for each mode. It is thus proportional to the one-third power of the 
frequency and to a constant value of which is now determined from 

the boiondary impedance. 

Introducing the asymptotic expression (25) into equation (21) to 
obtain the roots and writing = 1 + 2iT](o)/a) give 

.5 


^ . L ir \ _ 2 inBT^ 

i V t 2 + B(x2/n2) 


(27) 


For a pressure release, B = 0 and = jt[m - (l/^)J • This is the 
case which Pekeris treated and leads to Aq^ = I.85, A2 = 3 - 2 k, .... 

For a pc boundary B“^ = 1 + 2 i(o)/a) ^ 2 in and 


cot I& 


-|) = 


p 2 + 


(x2/2in5) 

This gives approximately 


.in/ 


1 + i(50)-2/5„-l/3e-«/6. , , 


V ' ”(■“ - - 


.in / 


for large values of n. Thus, in this case, the root depends slightly 
on the frequency but may always be ascertained for any given frequency 
by iteration, that is, by inserting the value [in - (l/^Q ®m on 

the right to obtain a second approximation, and so forth. 

For a hard boundary, B = 1 and equation (27) reduces to 
2 inT^ X 


cot 0 


+ (x2/n2) 


= - 2 inT^ 


= -6e 


i + 1 55/3(i)2/3e-i*/3. . . 
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This gives 0j_ = 0 .^ 69 , ^2 ~ 5-88^ ^-nd 6jj^ = n[^ - (5A^ } with m S’ 3; 

whence, A 2 ^ = 0.55 A 2 = 2.57* 

Numerical examples . - At large distances from the source the 
first mode is clearly the only surviving one. As an example, take 

a = 0.01 second”*^, which is a typical atmospheric value. At 5 OO cps 

the decay rates for the first mode are db/mile for Z = 0, 45 db/mile 

for Z = pc, and 12.5 db/mile for Z = c». As a second exainple, take 

a = 100 seconds"^, f = 10,000 cps, and Cq = 420 m/sec. The decay rates 

for this case are 8.6 db/foot for Z = 0, 8.4 db/foot for Z = pc, and 
2.6 db/foot for Z = 00 . 


Medium With a Constant Temperature Gradient 

In a medium with a constant temperature gradient the velocity func- 
tion is 

c(z) = Co[l - (pz/ToJ ^ 

where Tq is the absolute temperature at the groimd and p = dT/dz is 
the temperature gradient. To obtain a solution to actuation (l6) a method 
proposed by Langer (ref. 8) is adopted. Langer shows that 


^ + (q2 - 0)f = 0 

dz2 


where 



1 2 .( 9 ]^ 

2 Q 56wy 


has the exact solution 


F = ]J(u/q)Z2^^^(u) 


(28) 


in which 

u=r 


Q dz 


^1/5 


linear combination of Bessel functions 


first-order zero of Q 
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These solutions given in equation (28) are a good approximation to the 
actual solutions of equation (16) provided that q2 » 9 . Now, in this 

case, with = co^/ c(z) ^ the term 0 becomes (see appendix B; 

symbol 0 means "order of") 




+ 0 




lO) c H / 


with H = Tq/P. Also, it appears from equation (59) of the body of the 
text and equation (46) of appendix C that for high frequencies 




2/3 


The condition » 9 


may then be replaced by 


(coH/Cq) » 1 


(29) 


This restriction^ however^ is essentially the same as that already 
assumed in order for the wave eqiiation (eq. (l4)) to be valid, since 
f/a ^ ^^^/(jtCoP) where a = dc/dz . 

In terms of these functions, the solution given in equation (l9) 
takes the form 


^ ^ Ky pH, . r rpK 12 


( 2 ), 


( 2 ), 


(u/Q) 


Hi/5^^^u) 


dz 


where the eigenvalues are obtained from 






z=0 


(30) 


(31) 


Here, Ujj^(z = O) = u^, Ujj^(h) = u^^, and Ujj^(z) = u; similar expressions 
are obtained for Q. For the functions c(z) chosen. 
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(H = To/P) (52) 


and 


u 




tan-l(Q/k) - (H - z)Q 


(33) 


The substitution Q = k tan t reduces eqiiation (35) bo the form 


T 

COS T 


sin T + 



(3^+) 


with 


= 




-I3t 


Also, it is shown in appendix C that equation ( 51 ) may be approximated 

by 




Jl/3(vm) + J-l/3(^m) ^ -iZ ^ 


(35) 


For the higher modes the expression on the left is conveniently replaced 

by cot|vjjj - ()t/i^■^ . The determination of the roots of equation ( 31 ) 

then involves the simiiltaneous solution of equations (34) and (35)* 

Before solving these eq\jations, however, it is instructive to examine 
the general solution in the limit of high frequencies and so gain a 
better insight into the diffraction mechanism. 

Behavior of solution at high frequencies . - The quantities u, u^^, 

Q, and Qq appearing in the general solution, eqimtion (50), are related 
in equations ( 32 ) and (33) • From equation (34) it is clear that t is 
small in the high-frequency range. Accordingly, put 


= (a>/co)cos t = (a>/co)(l + bj^) 


( 56 ) 


Then, to the first order in 6 ^^, 
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Q « (cd/c) 


- 6„(c/oo)^(H/z)^/^ 


where the approximation holds provided (z/h) » Similarly, 


\ 1 + (z/H)(Co/c)2 


Introducing these expressions into equation (35) gives to the first order 


in b^ (note that tan"*^ (zc^/Ec^ 


_ -..-1/. Ml 1/21 


= sin (z/H) 

Uj^(z) » (a^H/co) |sin“^(z/H)^^^ - (z/h)^^^|i - (z/h^ 

^ j 

< 

sin"^(z/H)^^^ + (z/h)^^^|i - (z/h)| 


6m(cDH/co) 


The first term in this phase factor is recognized as the one obtained 


directly from ray acoustics, namely, o)|t'(z) - |^'(z)/c^ 

^ > 
sin“^(z/H) + (z/h)^^^[i - (z/H^ 


, in which 

•'k 


r’(z) = H> 


t'(z) = (2H/cQ)sin“^(z/H)^^^ 


(37) 


The quantities rQ = r'(h) + r'(z) and tQ = t'(h) + t'(z) are then 

the distance and travel time, respectively, of a ray from the source at 
height h to a point at height z on the shadow boundary. The small 
coefficient of the second term in the phase factor u is complex 

and leads to attenuation when r is greater than the shadow distance r^. 

This is seen by replacing the Hankel functions appearing in solution (30) 
by their asymptotic forms and by using equation (57) ^*or u and u^ to 

obtain the high-frequency approximation to the general solution 
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m=l 


ro^l/2^1o)|t-to- |(l+&m)(r-ro[j/co| 


(58) 


where 


2'[2 
^ " “2" 

F ‘ 




Here, 


= 


H 


0 


Hi/j'^>(u) 


-| 


dz 


2 Q-2g-2i[u-(5«/l2)] 


du 


The quantity !l^ is the pressure amplitude of the mth mode on the shadow 
homdary r = r^. In two dimensions the factor Tp/r is absent. 

The phase of expression (38) gives the time of arrival of each mode 
at a point in the shadow zone as 


r - r. 


t = to + 


+ 8m) 


This equation indicates that each mode travels along the downward leg of 
the limiting ray to the ground and then travels a distance r - rQ along 

the ground surface before finally reaching the point (r,z) by a path 
parallel to the upward leg of the limiting ray. This is just what would 
be expected from Fermat's principle, since the traversal of this path 
clearly takes the least time. The small frequency -dependent term 6^ 
leads to dispersion. Its imaginary part is responsible for the attenm- 
"tion which is seen to take place only over the horizontal stretch along 
the grovind surface. 

Rate of decay of sound in diffraction region .- The decay of sound 
in the shadow region due to diffraction alone is given essentially by 
the first mode as 8.68lmk3_ decibels per meter since the higher modes 
attenuate much more rapidly. The propagation constant k has been 
related to the temperature gradient p, the frequency f = oi/2«, and 
the impedance Z through equations {^k) and ( 35 )* Since b = Cq/cuH is 
small according to the criterion of equation ( 29 ), equation ( 34 ) reduces 
approximately to 

2 q -3 .j. 1 + bv = 0 

5 5 
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or, on solving for t, 


T 



(59) 


Either a high-impedance or a low-impedance approximation to equa- 
tion (55) te obtained by expanding the left-hand side in a Taylor’s 
series about the infinite or zero impedance values, respectively. Making 
use of recursion relations for the Bessel functions in obtaining the 
successive derivatives of f(v) and setting f(vo) = 0 then give, for 
the low-impedance boimdary. 


f(v) . - V - . 


bVr 


where Vq = 2.58 for the first root. Now, solve for v to obtain 
approximately 

.,2 


V « Vq - f - 


bv. 


(40) 


For a high-impedance boundary, put f(voo) = where Voo = O.685 
and define g(v) = l/f(v). Proceeding in the same way as before the 
equation 


g(v) 



is obtained; whence. 


V 


sL 

6Vo3 


(41) 


Finally, by expanding cos t, equation (56) may be replaced by 


Im k = - — — Im 

2 Cq 


(42) 


where Im k is the imaginary part of k. 

Low- impedance boundary . - An expression for t in terms of 
^ = Z/pqCq and Vq may be obtained by substituting from equation (59) 
into eqiiation (58) and writing equation (55) f = -i^r. This gives 
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T 


2 





Substituting this into equation {h2) gives 


Im k = 1.01 



1 

2 



( 45 ) 


This series holds in the range U P « (coTq/CoP) which under ordinary 

atmospheric conditions outdoors means |^| « 25 Tor frequencies of the 
order of 100 cps . By including more terms in the Taylor's expansion it 
can be shown that the coefficient of the next term in the series, pro- 

portional to Kl , vanishes. 

High-impedance boundary . - An expression for t in terms of 
T) = PqOo/^' and Voo may be obtained by substituting from equation (4l) 

into equation (59) and writing equation (55) as g = iq/r. This gives 



From equation (42) , 

Im k = ^ hi - 0 ) • • • (44) 

where 11 = ^iile-i0. This series holds in the range 

which imder ordinary conditions in the outdoor atmosphere means 
I ^ I » 25 for frequencies around 100 cps . 

In this connection it is of interest to note that the attenuation 
constant for a pres sure -release boundary (zero impedance) is 2.5 times 
that for a reflecting boundary (infinite impedance) . 


Typical numerical example . - Consider a point source located 10 feet 
above ground and driven at 5 OO cps in an atmosphere in which teir 5 )erature 
decreases with height at the rate of 0.5° C per meter. A receiver located 
at the same height passes into the diffraction region at a distance of 
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550 feet. Within this region the sound field hehaves like a damped 
cylindrical wave. In other words ^ apart from decreasing at the rate of 
5 decibels per doubling of distance because of divergence, the sound 
pressure level is observed, from egiiation (43), to diminish at the addi- 
tional rate of 78 decibels per 1,000 feet over a ground with a normal 
impedance of pc. If the impedance is lOpc this rate will be 58 decibels 
per 1,000 feet. From equation (44), a perfectly hard boundary results in 
a diminution of 55 decibels per 1,000 feet. Note that the major portion 
of the sound intensity decays exponentially with distance at a rate pro- 
portional to the two-thirds power of the ten 5 )erature gradient. 


LABORATORY MEASUREMENTS 
Propagation Chamber 


In the study of the sound field in the laboratory it is essential 
to obtain a high temperature gradient in the region of sound propagation 
in order to secure a shadow zone within a reasonable distance from the 
soimd source. To this end a three-dimensional propagation chsimber was 
first constructed. It consisted of a wooden box 8 by 4 by 3 feet with 
acoustically insulated side walls. In order to reduce the presence of 
air currents the temperature was made to increase with height. The 
’’ground” or hard boundary then consisted of a steel hot plate placed 
along the top of the chamber while water was circulated through coils 
in the bottom. The temperature at different heights within the chamber 
was measured by a series of mercury thermometers, which were shielded 
against radiation. Average temperature gradients of the order of 80^^ C per 
meter were set up within the chamber, but they were far from uniform, 
being greatest near the hot surface. Acoustic measurements were made by 
introducing a p\ire tone into the chamber at one end and then exploring 
the resultant field near the other with a condenser microphone (Altec 21B) 
with probe attachment. The probe consisted of a 4-foot Pyrex tube which 
was inserted into the box through holes in the top. The signal was 
amplified and measured with a nsirrow-band wave analyzer (General Radio) . 
Figure 5 is a plot of the intensity obtained as a function of depth, the 
labels ( 1 ), ( 2 ), and ( 3 ) referring to increasing distances from the 
speaker, which was driven at 8,700 cps. The presence of a shadow bound- 
ary is revealed in the rather sharp drop occurring in the intensity at 
a depth of about 10 inches. It is also noticed that the depth to this 
boundary increases with distance from the source as would be expected. 

The ’’shadow-zone distance” is in order-of -magnitude agreement with that 
predicted by ray theory (eq. (4)) for this case. 

The crude results obtained with this chamber pointed to the desira- 
bility of constructing an improved chamber in which a more stable and 
uniform temperature gradient could be set up. The final two-dimensional 
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propagation chaiuber, shown in figure 6, was found to be satisfactory 
in these respects. This chamber allows cylindrical waves to propagate 
between two plane parallel steel plates spaced a distance apart less 
than half a wavelength. Unwanted reflections are eliminated by means 
of a lining of sound -absorption material. The chamber is heated on the 
top and cooled at the bottom so that a vertical temperatiire gradient is 
maintained in the propagation region. 

The steel plates with the dimensions 8 feet by 5 feet by 5/16 inch 
are separated 5/8 inch, heat insulated on the outside, and supported 
vertically in a wooden frame. Five 600-watt heating units are placed 
along a 6-inch-wide plate which is bolted along the top of the chamber. 

The bottom is cooled by water circulating through a troiagh in which the 
plates rest . The teit5)eratiire is measured by means of merciory thermometers 
inserted into metal sleeves which fit into holes bored at 6-inch intervals 
in one of the steel plates. 

A cylindrical source of souind is simulated by a series of l/8-inch 
slots around the periphery of a. short length of pipe with a 3/^-inch 
bore and with a loudspeaker at one end. The pipe is inserted through 
the steel plates at a depth of 4 inches below the hot surface. The 
sound field is explored with a condenser microphftne with probe attach- 
ment. The probe consists of a 4-foot steel tube with a 3/16-inch bore 
and is inserted through holes at 10-inch intervals along the hot plate. 
These holes are mmibered from (l) to (9) starting from the source in 
figure 6, and plots made of measurements taken through them are nvimbered 
correspondingly in the graphs (figs. 9 , H. 15, and 14) . The microphone 
arrangement can be raised or lowered continuously by means of a motor- ^ 
driven elevator. Figure 7 is a photograph of the chamber and electronic 
instruments. Consideration was given to the question of what temperature 
distribution could be predicted in this experimental setup, allowance 
being made for the temperatirre dependence of the thermial conductivity of 
air. Since it is known that the conductivity A(T) increases with tem- 
peratin-e, it follows that the temperature gradient in the air space is 
nonmiform, having larger values in regions of low temperature. The 
actual distribution is given as the solution of the steady-state heat- 
flow equation 


subject to the appropriate boundary conditions . This solution is worked 
out in appendix D for a reasonable choice of A(t) , and the results are 
presented in figorre 8. The temperature midway between the metal side 
walls is plotted as a function of distance from the hot surface for dif- 
ferent values of the parameter 7 which is the ratio of Tq/ 3^ 

"height of the atmosphere" to the separation L of the side walls in 
which a constant gradient of temperature is assumed to exist. The curve 
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labeled 7 = 0 then corresponds to the case where the side walls have 
been removed to infinity. It is seen that in all cases the temperatiore 
gradient is essentially constant near the hot surface. 


Experimental Results 

An average temperature gradient of 210^ C per meter was maintained 
in the side walls during the sound measurements. Since it appeared that 
an inordinately long time would have to be allowed for the chamber to 
attain a uniform temperature gradient^ these measurements were usually 
taken after about 5 to 4 hours of heating, at which time the gradient 
at the top was still somewhat higher than the average one (about 250^ C 
per meter) . 

In the first series of measurements made in the new chamber the 
soiirce was driven continuously at various frequencies. The results of 
a typical set of measurements made at 9^250 cps are shown in figure 9- 
They are qualitatively similar to those obtained in the previous chamber. 
The sudden rise in the sound pressure level at the hot surface was always 
observed Just as the probe tip reached the hole orifice and is probably 
due to a diffraction effect. 

In all subsequent measurements the speaker was driven by pulses, and 
the resulting signal from the microphone, suitably amplified and filtered, 
was then displayed on an oscilloscope screen whose sweep was triggered by 
the pulse to the speaker. Relative pressure amplitude measurements were 
made by using an attenuator in series with the oscilloscope to maintain 
a constant pulse height on the screen. A block diagram of the electronic 
setup is given in figure 10. Figure 11 shows the results of a series of 
measurements taken in this way with the speaker driven by 8,000-cps milli- 
second pulses. The pulse length was clearly not short enough to eliminate 
interference caused by reflections from the hot surface in the normal 
region. 

To obtain a different picture of the distortion of the sound field 
by the temperature gradient, the relative times of arrival of the pulses 
at different points in the chamber were observed on the oscilloscope 
screen by noting the distance between the beginning of the sweep and the 
appearance of the pulse. Figures 12(a) and 12(b) show the results of 
such measurements taken in the chamber when cold and when heated, respec- 
tively. The numbers give the arrival times of the pulses in arbitrary 
units. Surfaces of constant phase have been drawn in by interpolation. 

The distortion of these surfaces in the heated chamber is apparent. 

In order to study the effect of an absorbing boundary on the sound 
field, a 2-inch layer of glass wool was placed directly under the hot 
surface. Figiore I5 gives the results of measurements taken with this 
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boundary using millisecond pulses at 8,000 cps, but with the chamber 
unheated. The increasingly steep drop in the pressure level as the 
boundary was approached at greater distances from the speaker suggests 
the shadow-zone effect. It may be explained q.ualitatively as caused 
by the increasing impedance mismatch between the wave and the boundary 
for decreasing angles of incidence and is predicted quantitatively in 
reference 2. 

The measurements reported in figure 14 were taken in order to illus- 
trate the additional sound attenuation brought about by the shadow-zone 
effect . The curves represent the differences in level obtained between a 
series of measurements taken at six different frequencies in the imheated 
chamber and the same series taken in the presence of a temperat\ire gra- 
dient. There is a noticeable increase in the slopes of these curves with 
frequency, altho\:igh for a given frequency the slopes remain fairly con- 
stant with distance from the source. A compaxison of these curves with 
those of figures 9 and 11 shows that the shadow zone has moved in closer 
to the source as would be expected since the introduction of a 2-inch 
layer of glass wool has reduced the source height. The curves afford a 
rather striking illustration of the fact that the influence of a tempera- 
ture gradient is quite negligible in the normal region but becomes pro- 
nounced in the diffraction region. 

The form of equation (58) indicates that a more direct comparison 
between theory and experiment wo\ild be possible if the pressiure measure- 
ments were plotted as a f\mction of distance from the source rather than 
depth, for such plots should yield straight lines in the diffraction 
region. Since the experimental setup was not adapted to taking such 
measirrements directly, they had to be obtained from the depth measure- 
ments by plotting successively points corresponding to the same depth. 

To justify this procedure, it was imperative to insure that the source 
output remained constant over the relatively long periods involved in 
moving the microphone and elevator from one access hole to the next. 

This was accomplished by introducing a monitoring microphone near the 
source. Its position is clear from figures 6, 7, and 10. When this was 
done, it was found possible to repeat all measurements within a variation 
of 1 decibel either way. 

In figures 15 to I7 the pressure level has been plotted as a function 
of distance from the source in meters. The curves in figure I5, repre- 
senting measurements taken at four different depths in the chamber, have 
been placed one above the other, and the approximate shadow bovmdary as 
suggested by the cxrrves themselves has been drawn in as a dashed line. 
Figures 16 and I7 are plots of pressure amplitude measurements in deci- 
bels taken at the source height with the hard and the absorbing bo\mdary, 
respectively, and at frequencies of 2,700, 5^000, 8,000, and 10,400 cps. 

In the shadow zone in both figures the level is observed to decrease 
linearly with distance from the shadow-zone boundary . 
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The theory given previously predicts that as the distance from the 
shadow boundary increases the behavior of the sound field is given by 
the first mode^ which decays at the rate of 8.68 Imkj_ decibels per unit 
of length. In order to compare the results in figures l6 and 17 with 
the theory^ assume for the temperatirre gradient the value prevailing 
at the top of the chamber^ namely, (3 = 250 ^ C per meter, since the 
rays are confined to this region. Correspondingly, = ^ 50 *^ C and 
Cq = 420 meters per second. The following decay rates are then obtained 
for the hard boundary and are compared with the measured slopes in 
figure l6: 


Frequency, cps 

Predicted, db/m 

Measured, db/m 

2,700 

8.8 

10 

5,000 

10.8 

14 

8,000 

12.7 

161 

10,400 

15.8 

19I 


If a pc impedance is assumed for the porous layer, the following values 
are obtained from figure 17: 


Frequency, cps 

Predicted, db/m 

Measured, db/m 

2,700 

19.5 

19 

5,000 

24.0 

25 

8,000 

28.1 

27 

10,400 

50.7 

29 


These measurements were repeatable to within ±1 decibel. The agreement 
is remarkably good for the soft -boundary case. Measured values with the 
hard boundary are somewhat greater than those predicted on the assumption 
of an infinite impedance, but this is probably to be expected since in 
practice the impedance is always lowered by the presence of a viscous 
boundary layer. 

Figure l8 is a logarithmic plot of these decay curves as a function 
of the frequency, which from the first terms of equations ( 43 ) and ( 44 ) 
should be straight lines of slope l/ 3 * This is again well confirmed in 
the soft-boundary case. The slope for the hard -boundary case is somewhat 
greater than l/ 3 - This is presimiably explained as a contribution from a 
finite admittance, due to viscosity and heat conduction, which is signif- 
icant for grazing angles of incidence. 
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CONCLUDING REMARKS 


A theoretical and experimental study of the sound field about a 
point source over a plane boundary in the presence of a vertical tem- 
perature gradient has been made. Results of this study indicate that 
the existence of a temperature gradient, which by itself leads only to 
a simple refraction of the sound rays, can in conjunction with the ground 
boundary result in the formation of a shadow zone. In the near field 
(normal region) the effect of the temperature stratification is slight. 

The field is essentially that due to a spherical source and its image in 
the boundary and, as such, depends strongly on the boundary impedance . 
However, the results of this study show that within the shadow zone or 
diffraction region the diminution of the sound intensity with distance 
is very much greater than that ca\)sed by ground and air absorption in a 
homogeneous atmosphere . In fact , in this region the sound behaves like 
a damped cylindrical wave, the intensity level of which, apart from 
cylindrical divergence, decays linearly with distance. The rate of decay 
has been obtained theoretically in terms of frequency, temperat\jre gra- 
dient, and ground ingjedance. The results, both theoretical and experi- 
mental, indicate that the major portion of the sound intensity in the 
shadow region decays exponentially with distance at a rate proportional 
to the one -third power of frequency and to the two -thirds power of tem- 
perature gradient . The effect of boundary impedance enters mainly 
thro\agh its resistive component. The rate of sound decay for a pressure- 
release or absorbing boundary (zero impedance) is found to be 2.5 times 
that for a reflecting boundary (infinite impedance) and for low impedances 
is largely independent of the reactive component. 


Massachusetts Institute of Technology, 

Cambridge, Mass., December 12, 195^. 
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APPENDIX A 

APPROXIMATION OF A SLOWLY VARYING MEDIUM 


The general linear wave equation is (ref, 8) 

^ = cV 0 + (Vrto)V0 

.../s 

is assiimed^ where 7 is the ratio of specific heats and r is a con- 
stant, the following equation may he obtained: 


It 


o 



p(7-l) 


Whence, for a stratified medium. 


Vrto = - - 

° dz o p p 


Now, = ZE; hence, 
P ’ ^ 


- 1 <=(? 


i)£i 


Thus, 


V«o 


2cc ' 


If it is assumed that 0 « ^ it follows directly that the 

condition 


c^20 » (Vjto)V0 


is equivalent to 
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APPROXIMATION TO WAVE EQUATION 


To evaliiate 


where 


compute first 


0 = 


5fQ'\^ 1 Q" 

k\Q/ 2 Q 

5 /Q' 
56 Vu. 

a" = 2^ - 


u = r Q d2 
^0 

r 

c = Co(l - z/H)^/^ 

c' = -Cq^/pcH 

0 

-=t 

0 

0 

1 

II 

"0 

r 2 

2 2 

At 03 Co'=^ 

^ *" — E — 
2c^QH 


4 E 

Q" = + 



u = -(c/Co)^QH + 


2^ti . ^Q^H/cQ%jtan"^(Q/k) 


When Q/k is small, 
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U— > 


c^QH 




_k 2 + ^ _ + 

c2 3Vck/ 5\ck2; 


Substitute = — (l - 7 ) where y = Then, 


u = + i r ± ,2 * . . . 

5Cq 0 )' 


2 -2 \ 5 ■ 10 


and 




By combining these, the following equation is obtained: 


e = -(5/16) (co/c)Vh 2 + o(qc//(dch)' 


+ • 
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appendix C 


SIMPLIFICATION OF EIGENVALUE EQUATION 


In order to find the eigenvalues the roots of equation (51) 

must first be determined. The right-hand side of this equation is 


d^m d 
dz dz 


/u\l/2^ (2) ^ 

fe) 


6W " 2^Q^ Q 


z=0 




so that 


iCDp 


1 ^ 
F dz 


z=0 


Qo 1 V 

2 Qo 


6 u^ 


+ Q^ 




H 




(45) 


where the primes denote differentiation with respect to z. As was shown 
in appendix B_, the values of Q/u given in equations (32) and (33) 
be expanded in terms of the small quantity 7^ as 


%/^o = (^ ■ 5 '^O ■ ^ ' • •) 

where 

= (^o%^Y = (3coUo/anH)^/5 

Flurthermore, 
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from which it follows that 




dz 


/ z=0 


i H + 




Ti ( 2 ),,^ - 

%/5 (^) 


Since the first term on the right side is very small compared with the 
second, it is justifiable to approximate equation (it-5) 


^ -IZSo „ ^- 1 ,\ (46) 

“®° mo/ 


The last step follows from the relations 



H ^^^(z) 

V 

= (i/sin vn)[j_^(z) - J^(z)e^'^ 






p 

f o 1 p\ 

ill/2 


Now, Qq = 

) - 

k is itself a function of 

through 

equation (55) 

taken at 

z = 0. With the substitution Qq = k tan t it 

follows that 






Qo = (“Ao)sin T 

(47) 


which together with equation (46) gives equation (35) • 
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APPENDIX D 

TEMPERATURE DISTRIBUTION IN A RECTANGULAR ENCLOSURE 
The steady-state heat -flow equation is 


div[x(T) grad = 0 


(48) 


where X(t) is the conductivity and T, the temperatiire . This equation 
is nonlinear and is not easily solved directly. Storm (ref. 9) shown 

in a study of heat conduction in metals that in the one -dimens ional case 
it can be linearized provided that 




dT 


log 


r 


pe- 


ls constajit, where p is the density and the specific heat at con- 

stant pressure. However, this condition does not seem to be satisfied 
in the case of a gas for any reasonable choice of A(t) . 


Kinetic theory based on a rigid-sphere model predicts that for 
molecules repelling each other as the inverse vth power of the distance 
the coefficients of viscosity and heat conduction should vary as T®, 

where s = i + — - — (ref. lO) . This tj'pe of dependence is in good 

2 V - 1 


agreement with experimental results obtained by Stops (ref. ll) who 
measured the thermal conductivity of air over a wide temperature range 


(0° to 1,000° C). Accordingly, put A = ^o|^y 
to obtain 


in the equation (48) 


V^T + ^(grad T)^ = 0 


(49) 


It is recognized that this equation can be linearized by the transformation 


T = t' 


1 

s+1 


(50) 
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so that 


v2t = 0 


(51) 


In the experimental setup the high temperature is on the top of the 
chamber in order to diminish the effect of convection of air which is 
disregarded here. Furthermore, radiation is neglected by assuming that 

oT^ « grad T (a, Stefan-Boltzmann constant; conductivity of 

the metal) . Moreover, since the ratio of the conductivity of metal to 
air is on the order of 5,000 for steel, the heat transfer between the 
air and the side walls will have a negligible effect on the temperature 
distribution. To a first approximation, the temperature gradient within 
the side walls can therefore be assiomed to be constant. The problem is 
then to find a solution to equation (^9) satisfying the boundary condi- 
tions T(x,0) = Tq and = "^o “ variables are changed 

as 

Tl/(=+l)(x,y) = 

'^O 

the problem reduces to that of obtaining a solution to Laplace's equation 
in a bounded two-dimensional region, 0 = y = '^o/P^ 

boundary conditions now being t(x,0) = 1, t(x,Tq/P) = 0, and 

T(t|,y) = (|Sy/To)"-^l 

Set r(x,y) eq\ial to a linear combination of fvmctions 
Tp(x,y) + T2(x,y), where and T2 are separate solutions to 

Laplace's equation satisfying the boundary conditions, tq^(x,0) = 1, 

■'■l(^>'^o/P) = T3^(-L/2,y) = 0 and T2 (x, 0) = T2 (x,Tq/P) = 0, 

T2(tL/2,y) = (Py/To)^'*’^ Expressing the solution as a Fo\rrier series 

in the usiial way gives 


"^1 


4"^ (-l)”^ sinh(2m + 1)«7 Ti /„ . n \ / h\ 

- > 1 ) . ' f cos(2m + l)(rtx/L) 

Jt m=0 2m + 1 sinh(2m + Ijity 


(52) 
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^2 


00 



m=0 



s+1 


X 


cosh(rmtpx/To) 

sin imtx dx sin imtT) 

cosh(imt/27) 


where 


(53) 


3L 


7 



- n 


Use of the identities^ 

ey - e~y ^ Q-y+y_Q-y slnh(7 - y) 
e7 _ Q-y ~ Sinh 7 


and 


i = i tan-1 

m=0 ' 


where Re indicates "real part of," improves the convergence of 
which then becomes 


T-| = — |— tan"^ 
1 n 2 


cos rtx 


Ik. 


sinh rt7(l - t]) 


m=0 

At X = 0, 


( sinh(2m -H l)jt7(l - t]) 
2m -f 1 sinh(2m + l)rt7 


cos(2m + 1 ) (jtx/l) 


Ti(0,T]) = i 


H 

m=0 


tan-le-«’'(l-^^ - 

sinh(2m -1- l)jt7(l - t]) 


sinh(2m -1- l)«7 


2m + 1 
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In the case of T2 it is readily verified that the value of the 

integral which appears in the Fourier coefficients does not depend crit- 
ically on the value of s for n > 2 provided 0 S s ^1. Accordingly^ 

s = 0 may he taken for n > 2 and the integral may he evaluated hy 

Simpson’s rule for n = 1 and n = 2 with s = 0.75* Then, for T2 

at X = 0, the following equation is obtained: 


"2 =! 


0.656 


sin itT] 
cosh(jt/27) 


- 0.512 


sin 2 jit] 


(- 1 ) 


ra +1 


sin iMT] 


cosh(rt/7) m =5 m cosh(imt/27) 


Figure 8 shows curves of T = (T]^ + for different values 

of the parameter 7 = Tq/PL. The curve labeled 7=0 corresponds to 
the case in which the side walls have been removed to infinity and is 
equivalent to a situation in which the side walls are replaced by 
nonconductors . 
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Figure 1.- Illustration of mechanism of shadow-zone formation. 
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Figure 2.- Distance Tq to shadow-zone boundary at source height h in 
temperature gradient p. Temperature at ground is Tq (°K). 
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Figure 5-- Ray diagram illustrating computation of intensity distribution 
around a source in a stratified medium by considering geometrical 
divergence of rays. 
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Figure 4 .- Computed intensity I from a spherical and a dipole source 
in temperature gradient p referred to intensity Iq from same 
sources in homogeneous medium plotted versus distance parameter pr/To- 
Tq is absolute temperatiire at source elevation. 
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Figure 5-“ Intensity level in three-dimensional chamber plotted versus 

depth from hot surface at distances of (l) 6 feet, (2) 6 feet 6 inches, 
and (5) 7 feet from speaker. 
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Figure 6.- Schematic views of two-dimensional propagation chamber. 








Figure 7 .- Photograph of two-dimensional propagation 

equipment . 
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Figure 8.- Absolute temperature ratio t/Tq in air midway between two 
vertical side walls in which a constant temperature gradient p 
exists plotted versus height parameter py/T© for different separa- 
tions L = To/py of side walls. 
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Figure 9-- Measured sound pressure level as a function of depth from a 
hard boundary taken in chamber in presence of a temperature gradient 
at ( 7 ) 66 inches^ ( 8 ) 76 inches^ and ( 9 ) 86 inches from speaker which 
is driven by 9 ^ 250 -cps millisecond pulses. 
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Figure 10.- Block diagram of electronic setup. 
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Figure 11.- Measured sound pressure level as a function of depth from a 
hard boundary taken in chamber at (l) L inches, (4) 35 inches, 

(6) 55 inches, and (8) 76 inches from speaker, which was driven 
by 8,000-cps millisecond pulses. Temperature gradient wsls about 
210° C per meter. 
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(a) Without a temperature gradient. 
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(b) With a temperature gradient of 210° C per meter. 

Figure 12.- Surfaces of constant phase in two-dimensional propagation 

chamber . 
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Figure 15.- Measured sound pressure level as function of depth from 
absorbing boundary taken in chamber at (l) 4 inches, (4) 55 inches, 
(6) 55 inches, and (8) 76 inches from speaker which was driven by 
8,000-cps millisecond pulses. Temperature gradient was about 210° C 
per meter. 
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Figure l4.- Additional attenuation brought about by presence of temperature 
gradient plotted for different frequencies as a function of depth. 

Niombers on curves correspond to distances from speaker shown in figure 7* 
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Figure 15-- Measured sound pressinre levels at various depths in chamber 
as a function of horizontal distance from speaker driven by 9,000-cps 
millisecond pulses over absorbing boundary with vertical temperature 
gradient. Dashed curve indicates shadow-zone boundary. 
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Figure l 6 .- Measured sound pressure level as a function of horizontal 
distance from source over a hard boundary vith a vertical temperature 
gradient. 
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Figure I7.- Measured sound pressure level as a function of horizontal 

distance from source over absorbing boundary with a vertical temperature 
gradient. 
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Figure l8.- Measured attenuation in shadow zone as function of frequency, 
(a) indicates soft boundary; (b) indicates hard boundary. 
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